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INTRODUCTION 

STABL  is  a  computer  program  for  the  general  analysis  of  slope  sta- 
bility by  a  two  dimensional  (2-D)  limiting  equilibrium  method.  The  pro- 
gram uses  the  method  of  slices  to  analyze  the  slope  and  calculates  the 
factor  of  safety  (FOS)  according  to  the  simplified  Janbu  or  Bishop 
methods  for  non-circular  and  circular  surfaces,  respectively.  A  unique 
feature  of  the  program  allows  random  surfaces  to  be  generated,  allowing 
the  user  to  determine  the  critical,  minimum  FOS  more  easily. 

DEVELOPMENT  OF  STABL 

The  2-D  computer  program  STABL  was  developed  at  a  time  when  most 
highway  agencies  analyzed  slope  stability  using  two  common  techniques: 
computer-aided,  grid-type  circular  searches;  and  block  analyses  for  sim- 
ple and  specified  surfaces.  Circles  were  often  assumed  to  be  the 
appropriate  shape  for  potential  failure  surfaces  simply  because  there 
was  no  other  shape  which  could  be  used  for  computerized  searching. 

In  the  last  decade,  improvements  in  2-D  slope  stability  analysis 
have  proceeded  in  several  directions;  one  of  these  is  contained  in 
STABL,  in  the  form  of  computerized  searching  with  non-circular  shapes. 
The  non-circular  routines  RANDOM  and  BLOCK  were  first  reported  by  Siegel 
(1975a)  as  well  as  a  random  (as  opposed  to  a  grid)  type  search  with  cir- 
cles (CIRCL2).  Favorable  comparisons  of  the  FOS  values  generated  by 
STABL  with  those  for  the  same  surfaces  by  other  methods  of  slices  were 
reported  by  Boutrup  (1977). 


STABL  was  placed  on  line  for  routine  use  in  1976  by  the  Indiana 
Department  of  Highways  (IDOH),  and  after  being  reported  in  the  open 
literature  (References  11,  2,  3,  4,  12),  the  program  began  to  be  adopted 
by  many  agencies.  STABL  has  been  modified  in  minor  ways  over  the  past 
eight  years,  and  users  of  the  program  have  helped  greatly  in  debugging 
operations.  The  present  version  of  STABL  is  called  STABLA. 

CAPABILITIES  OF  STABL 

STABL  is  the  only  known  program  to  contain  searching  routines  for 
shapes  other  than  circles.  Circular  grid  search  routines  are  common; 
STABL  searches  by  generating  circles,  with  randomly  selected  radii  and 
centers  limited  by  user-defined  parameters  as  explained  later.  A  number 
of  programs  can  determine  the  FOS  for  any  specified  non-circular  shape, 
but  only  STABL  can  both  generate  general  shapes  and  determine  their  FOS 
values. 

The  2-D  program  uses  a  simplified  method  of  slices  and  thereby 
minimizes   iterative  procedures  in  the  solutions.   The  single  limitation 

on  boundary  geometry  is  that  there  be  no  vertical  or   overturned  (>90  ) 

surfaces.   Subsurface  boundaries  may  demonstrate  any  degree  of  natural 

complexity,  and  up  to  ten  piezometric  surfaces  may  be  specified.  Also, 

options   allow  boundary  and  pseudo-static  earthquake  loadings  to  be  con- 
sidered in  the  analysis. 

The  circular  potential  sliding  surface  is  appropriate  when  the  sub- 
surface materials  are  grossly  homogeneous  and  isotropic.  Block-type  sur- 
faces are  probable  when  weak  strata  are  present,  and  the  critical  sur- 
faces  tend  to  have  a  maximum  length  within  these  strata.   Thus,  we  have 


cited  two  cases  in  which  simplifying  assumptions  are  appropriate.  In 
the  first  case  it  is  proper  to  assume  a  circular  shape  and  iterate  for 
the  position  of  the  critical  surface.  In  the  second  case,  the  positions 
of  the  major  portion  of  the  sliding  surface  may  be  assumed,  but  the 
entry  into  the  weak  layer,  the  length  of  sliding  surface  in  this  layer, 
and  the  departure  from  it  need  to  be  iterated. 

A  far  more  common  occurrence  is  the  profile  where  neither  shape  nor 
position  may  be  safely  assumed  and  the  analysis  should  both  generate  and 
compute  FOS  values  for  a  wide  variety  of  shapes  and  positions.  The 
non-circular  surface  generating  routine  RANDOM  is  well  suited  to  this 
requirement. 

In  using  the  searching  routines  CIRCL2,  BLOCK  (BL0CK2),  and  RANDOM, 
the  positions  of  the  surfaces  are  not  obvious.  With  RANDOM,  the  shape 
generated  is  a  further  uncertainty,  and  plotting  routines  which  resolve 
both  uncertainties  are  necessary.  These  routines  can  both  show  the  sub- 
surface space  searched  and  the  portion  of  that  space  occupied  by  the 
more  critical  surfaces. 

COMPUTER  REQUIREMENTS 

STABL  is  written  in  FORTRAN  IV  and  has  been  successfully  used  on 
the  following  computer  systems: 

1.  CDC  6500-6600 

2.  IBM  360/370 

3.  VAX  11/780 


With  some  modifications,  STABL  may  also  be  converted  for  use  on  a 

micro-computer   (Ref.   10).  However,  the  versions  available  from  Purdue 

University  should  be  compatible  with  most  systems  which  can  compile  the 
FORTRAN  IV  language. 

On  the  CDC  system,  the  program  requires  a  field  length  of  approxi- 
mately 70,000  to  80,000  words  for  storage  and  execution.  The  CPU  time 
for  "average"  runs  is  generally  16  to  32  seconds  for  approximately  100 
circular  surfaces  and  may  be  up  to  128  seconds  for  100  irregular  sur- 
faces. However,  the  cost  of  each  run  is  directly  proportonal  to  the 
complexity  of  the  surface  and  sub-surface  geometries  and  will  be  a 
unique  feature  for  each  individual  analysis. 

DESCRIPTION  OF  STABL 
Assumptions 

STABL  assumes  that  the  instability  to  be  prevented  would  be  two- 
dimensional.  In  reality,  all  sliding  failures  must  be  3-D,  with  the 
end/edge  resistance  furnishing  additional  safety  against  instability. 
For  more  quantitative  information  on  the  comparison  of  (F0S)__  to 
(F0S)2D,  see  Chen  (1981)  and  Lovell  (1982).  In  general,  (F0S)3D  > 
(F0S)~_,  but  the  difference  may  be  small,  and  in  certain  special  cases 
(F0S)~  >  (F0S)3  .  where  the  stability  problem  is  perceived  to  be 
definitely  3-D,  the  engineer  is  encouraged  to  use  the  BL0CK3  or  LEMIX 
codes  of  Chen  (1981). 

STABL  uses  simplified  methods  of  slices  for  determination  of  F0S. 
The  alternative  requires  solutions  with  extensive  iteration  and  the  con- 
sequent problems  of  nonconvergence  in  these  iterations.   Boutrup   (1977) 


has   shown  that  the  simplified  methods  after  Janbu  and  Bishop  give  rea- 
sonably precise  values  of  FOS. 

Agencies  staffed  with  appropriate  mathematical  and  software  skills 
can  insert  any  desired  slices  solution  into  the  program. . .simplified  or 
total  equilibrium.  STABL  is  a  stability  analysis  system,  of  which  the 
method  of  slices  detail,  is  a  small  part. 

The  selection  of  a  center  of  moments  for  the  slices  analysis  is  an 
intriguing  point.  In  the  simplified  approaches,  the  free  body  is  not 
iterated  into  equilibrium,  and  accordingly,  the  FOS  value  is  peculiar  to 
the  center  selected.  This  is  true  even  for  the  circle,  where  the  circle 
center  is  arbitrarily  selected  in  the  simplified  Bishop  method.  For 
other  shapes,  there  is  usually  no  "center"  to  select  for  moments.  After 
much  study  of  this  question  (Carter,  1971;  Siegel,  1975a,  Boutrup, 
1977),  the  circle  center  is  used  for  CIRCL2,  and  a  very  long  moment  arm 
is  used  for  BLOCK,  BL0CK2,  and  RANDOM.  The  latter  choice  means  that 
these  noncircular  surfaces  are  analyzed  with  the  same  slice  assumptions 
as  the  simplified  Janbu  method. 

Both  of  these  slices  methods  contain  the  FOS  value  in  implicit 
form,  which  requires  an  initial  iteration  for  FOS.  This  iteration  is 
greatly  expedited  if  it  is  entered  with  a  reasonable  value  of  FOS.  This 
can  ordinarily  be  accomplished  by  means  of  a  Taylor  type  chart  (Refer- 
ence 13),  where  the  real  problem  has  been  suitably  simplified  to  permit 
approximate  chart  solutions. 

STABL  values  may  be  checked  for  a  specific  failure  surface  in 
several  ways.   CIRCL2   should  yield  about  the  same  FOS  (for  the  same 


circle)  as  any  other  computerized  analysis  for  circles.  To  determine 
that  this  is  indeed  the  case,  the  new  user  of  STABL  can  run  CIRCL2  in 
parallel  with  his  present  method.  BLOCK  or  BL0CK2  can  be  checked 
approximately  (for  a  specific  block)  either  manually  or  perhaps  by 
existing  charts.  RANDOM  is  amenable  to  approximate  manual  checks. 

Input  Options 

Numerous  options  are  available  for  use  with  the  STABL  program 
according  to  the  user's  specific  requirements.  These  options  are 
defined  by  command  words  (e.g.,  PROFIL,  EQUAKE,  etc.)  for  easy  recogni- 
tion, especially  when  checking  input  data.  Different  commands  activate 
or  deactivate  different  portions  of  the  program,  which  also  allows  the 
user  to  control  and  perform  the  analysis  economicaly.  Generally,  these 
commands  may  be  categorized  into  four  major  divisions: 

1.  Surface  and  subsurface  geometry,  with  command  words 

PROFIL 
LIMITS 

2.  Subsurface  profile  parameters,  with  command  words 

SOIL 

ANISO 

WATER 

3.  Boundary  loads,  with  command  words 

LOADS 

TIES  (Available  in  STABL4  only) 

EQUAKE 

4.  Analytical  methods,  with  command  words 

SURFAC/SURBIS 


CIRCL2 
RANDOM 
BLOCK 
BLOCK2 

The  command  PROFIL  is  used  to  label  input  data  of  the  slope 
geometry  and  any  subsurface  boundaries  forming  the  soil  strata.  Up  to 
one  hundred  different  boundaries  can  be  specified.  If  the  user  is  try- 
ing to  confine  critical  surfaces  within  certain  zones,  the  command  LIM- 
ITS may  be  used.  This  prevents  the  generation  of  critical  planes  beyond 
the  defined  limiting  boundary,  such  as  competent  rock  underlying  the 
slope.  Of  course,  one  could  simply  introduce  the  rock  as  a  soil  type 
with  a  high  strength,  which  would  effectively  lead  to  higher  factors  of 
safety  (FOS)  for  surfaces  passing  through  the  rock.  However,  the  LIMITS 
option,  by  eliminating  the  surfaces  which  would  have  passed  through  the 
rock,  is  more  efficient  since  computer  time  is  not  wasted  in  computing  a 
high  FOS.  Additionally,  an  upper  limiting  boundary,  which  will  force 
planes  downwards,  may  also  be  specified.  This  is  useful  for  stability 
analysis  of  retaining  structures  such  as  sheet  piles,  anchored  bulk- 
heads, reinforced  earth  walls,  etc. 

In  order  to  define  subsurface  conditions,  the  command  SOIL  is  used 
to  assign  for  up  to  20  soil  types  the  soil  parameters  which  include  unit 
weights,  c-<{>  parameters  and  pore  pressures.  If  the  user  feels  that  the 
strength  properties  are  significantly  anisotropic,  the  command  ANISO 
allows  one  to  define  the  strength  parameters  in  up  to  10  sectors,  the 
appropriate  values  being  selected  according  to  the  plane  of  interest. 
The  WATER  option  is  used  to  define  up  to  10  different  piezometric  levels 


within  the  slope  geometry.  These  can  be  successfully  used  to  assign 
perched  water  tables  and  also  to  simulate  pore  pressure  distributions 
within  the  soil  strata. 

External  boundary  loads  which  may  be  applied  to  the  slope  are 
defined  by  the  command,  LOADS.  This  allows  input  of  up  to  10  loads,  as 
a  surcharge  intensity,  in  any  direction.  A  boundary  load  of  varying 
intensity  may  be  approximated  by  a  set  of  equivalent  uniformly  distri- 
buted loads  abutting  each  other.  Tieback  or  concentrated  loads  applied 
to  the  slope  are  defined  by  the  command,  TIES.  This  option  allows  input 
of  up  to  10  horizontal  or  inclined  concentrated  loads  applied  to  the 
profile  surface.  The  TIES  option  is  fully  described  in  Appendix  F  and 
is  available  only  in  the  most  recent  version  of  STABL,  STABL4.  For  a 
pseudo-static  earthquake  analysis,  the  command  EQUAKE  is  used  to  assign 
vertical  and  horizontal  coefficients  to  simulate  the  design  earthquake. 

Four  different  methods  of  analyzing  the  stability  of  a  slope  are 
available  for  use  in  STABL.  The  commands  SURFAC  (for  Janbu  Analysis)  or 
SURBIS  (for  Bishop  analysis)  are  used  for  determining  the  FOS  for  a 
specific  critical  failure  surface  defined  by  coordinates  input  by  the 
user.  However,  if  the  user  is  at  an  initital  investigation  phase,  the 
commands  RANDOM  (irregular  surfaces)  and  CIRCL2  (circular  surfaces)  may 
be  specified  to  utilize  the  program's  unique  capability  to  generate 
critical  failure  surfaces  in  a  random  manner.  The  user  only  needs  to 
define  the  initiation  and  termination  limits  of  the  failure  surfaces. 
If  block-type  analysis  is  also  required,  the  BLOCK  or  BL0CK2  options 
allow  the  user  to  generate   random  blocks  for  analysis.   For   these 


blocks,  the  user  can  manipulate  the  size  of  the  block  zones  to  arrive  at 
the  minimum  FOS. 

Thus,  by  specifying  the  geometry  and  soils  of  the  slope,  the  user 
can  perform  numerous  operations,  using  three  different  types  of  critical 
failure  surfaces,  to  arrive  at  the  most  probable  value  for  the  FOS  for 
the  slope.  Additionally,  all  the  commands,  except  for  PROFIL  and  the 
analysis  commands,  may  be  turned  "off"  or  "on"  which  allows  use  of  the 
same  input  data  for  the  entire  stability  analysis,  without  having  to 
create  multiple  data  sets. 

Data  is  input  as  "free-format",  beginning  in  the  first  column  (left 
justified),  with  items  being  separated  by  one  space.  This  facilitates 
data  entry  as  it  does  not  restrict  the  user  to  certain  column  widths 
which  sometimes  lead  to  format  errors.  Also,  STABL  is  compatible  with 
any  unit  systems  which  the  user  may  prefer  to  utilize  for  analysis.  Pro- 
viding consistent  units  of  length  and  mass  are  used,  the  program  will 
function  without  error.   For  example  one  may  use  the  following: 


Coordinates: 
Unit  Weights: 

Strength  Intercept: 

Pore  Pressure: 


feet  (ft) 

pounds /cubic  foot 
(pcf) 

pounds/square  foot 
(psf) 

pounds/square  foot 
(psf) 


metres    (m) 

kilograms/cubic  metre 
(kg/in   ) 

kilograms  /square  metre 
(kg/in   ) 

kilograms/square  metre 
(kg/ni   ) 


Unit  Weight    of  Water:        62.4   pcf 


1000   kg/m" 


10 
Searching  Routines 

Each  slope  to  be  analyzed  is  made  left-facing  and  placed  in  a  first 
quadrant  position  so  that  all  coordinates  will  have  positive  values.  In 
order  to  make  the  computer  generation  and  search  reasonably  efficient,  a 
number  of  restrictions  are  placed  on  the  surfaces: 

(1)  They  should  enter  and  exit  the  surface  boundary  at  reasonable 
positions.  This  is  accomplished  by  assigning  initiation  and  termination 
limits.   See  Figure  1. 

(2)  They  should  not  go  deeper  than  a  depth  limit,  usually 
represented  by  a  hard  layer.   See  Figure  1. 

(3)  They  should  be  composed  of  straight  line  segments,  not  so  long 
as  to  fail  to  represent  real  changes  in  surface  shape,  and  not  so  short 
as  to  produce  "kinkyness"  in  shape.  Guidelines  have  been  worked  out 
through  experience. 

(4)  The  generated  surfaces  shall  appear  to  be  kinematically  accept- 
able, i.e.,  when  examined,  movement  along  them  appears  logical. 

(5)  While  surfaces  are  generated  in  a  random  fashion,  there  must  be 
some  bias  inserted  to  obtain  acceptable  shapes. 

Figure  2  shows  schematically  how  the  first  line  segment  is  gen- 
erated. Without  appropriate  upper  and  lower  limits,  the  surface  could 
immediately  assume  an  illogical  shape  and  position.  The  direction  of 
the  first  line  segment  is  determined  within  the  limits  by  a  random 
number  generator  in  the  program.  As  shown  in  Figure  3,  the  second  seg- 
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FIGURE  2.     GENERATION  OF  FIRST  LINE   SEGMENT 
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ment  of  the  potential  failure  surface  has  an  assigned  upper  limit  and  a 
lower  limit  which  is  generated  in  a  random  but  somewhat  biased  direc- 
tion. The  direction  of  the  segment  is  again  randomly  generated  within 
the  limits. 

If  a  circle  is  being  generated  (CIRCL2),  it  has  been  defined  by  the 
random  generation  of  the  first  two  segments  (chords).  The  program  will 
complete  the  generation  of  the  circle  and  determine  its  FOS.  For  the 
irregular  surface  generation  (RANDOM),  random  choices  continue  to  be 
made  as  for  the  second  segment  until  the  surface  is  totally  generated; 
it  is  then  analyzed  for  the  FOS. 

Some  surfaces  will  violate  the  depth  or  termination  limits  (Figure 
1)  and  will  be  discarded.  Surface  generation  is  accomplished  in  a  batch 
mode,  with  the  generation  continuing  until  the  desired  number  is  accom- 
plished, or  until  a  large  number  of  attempts  have  been  made  without 
total  success . 

The  searching  procedures  described  above  are  appropriate  for  CIRCL2 
and  the  RANDOM,  but  may  not  be  economical  where  weak  layers  govern  the 
location  of  critical  sliding  surfaces.  In  such  cases,  grided  "boxes" 
are  superimposed  on  the  weak  layers.  As  shown  in  Figure  4,  boxes  are 
specified  to  overlie  a  weak  layer  wherever  that  layer  changes  in  direc- 
tion. 

Points  are  randomly  chosen  within  the  specified  boxes,  and  con- 
nected by  straight  lines.  The  surfaces  are  completed  at  the  ends  by 
random  generation  (BLOCK)  or  by  assuming  that  these  ends  are  simple 
Rankine  type  surfaces   (45 ±  y),   as   in  BL0CK2 .  The  routine  BLOCK  is 
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Weak  Layer 


Search  in  Irregular  Weak  Layer 


FIGURE    4.     SLIDING  BLOCK  GENERATOR  USING  MORE  THAN  TWO  BOXES 
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illustrated  in  Figure  5,  where  the  random  limits  cause  the  active  side 
segments  to  be  >-  45  ,  and  the  passive  side  segments  to  be  <-45  . 
BL0CK2  is  somewhat  simpler  than  BLOCK,  and  has  tended  to  be  used  more. 
Since  the  boxes  can  degenerate  to  points  or  lines,  the  program  can  be 
used  for  cracks,  joints,  or  fissures. 

As  emphasized  further  in  the  next  section,  interaction  between 
searching  routines  and  plotting  routines  is  required.  When  arbitrary 
assumptions  are  made  with  respect  to  selection  of  a  searching  mode,  ini- 
tiation zone,  termination  zone,  depth  limit,  box  position  and  the  like, 
the  plotting  routines  have  a  capability  to  reveal  this. 

Plotting  Routines 

The  use  of  plotting  routines  is  essential  to  the  random  generation 
of  surfaces  in  STABL.  These  plots  will  indicate:  -  the  shapes  generated; 
the  subsurface  space  searched;  and  the  positions  of  generated  surfaces 
which  have  the  lower  values  of  FOS.  Surfaces  are  generally  generated  in 
batches. 

The  plot  of  all  surfaces  in  a  batch  which  involved  10  irregular 
surfaces  from  each  of  5  initiation  points  is  shown  in  Figure  A2  in 
Appendix  A.  Also,  an  accompanying  plot,  which  depicts  ten  critical  sur- 
faces with  the  lowest  FOS,  is  shown  in  Figure  A3.  Depending  upon  the 
practical  problem  involved,  the  engineer  may  decide  to  change  the  ini- 
tiation limits,  terminiation  limits,  or  depth  limits  to  "explore"  other 
areas  for  a  more  critical  surface.  He  may  even  decide  to  use  different 
generation  modes,  i.e.,  BLOCK  (BL0CK2)  or  CIRCL2. 
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In  any  reasonably  complex  problem,  multiple  man-machine  interac- 
tions are  required,  and  are  anticipated  in  the  procedures  of  STABL.  The 
engineer  should  avoid  promising  the  final  results  of  an  analysis  based 
upon  the  computer  "turn  around"  time  for  a  single  run. 

Error  Messages 

In  a  comprehensive  program  such  as  STABL  it  is  likely  that  errors 
will  be  committed  in  data/analysis  input  and  commands.  Most  data  errors 
will  be  detected  by  the  program  and  the  user  is  notified  by  an  error 
message  and  a  print-out  of  the  line  of  erroneous  data.  These  error  mes- 
sages are  documented  in  the  User  Manual  (Ref.  10),  and  allow  the  user  to 
quickly  correct  mistakes  with  a  minimum  amount  of  time  spent  searching 
the  input  data.  Additionally,  data  are  checked  and  verified  continu- 
ously during  execution  of  the  program  to  minimize  the  risk  of  incorrect 
or  misleading  output. 

Recent  Changes 

Appendices  E  and  F  contain  modifications  effected  in  STABL  during 
the  past  year.  One  of  these  is  a  change  in  the  pseudo-static  earthquake 
option,  which  eliminates  the  increase  of  pore  pressures  for  an  effective 
stress  analysis.  Appendix  F  contains  an  improvement  which  makes  it  more 
convenient  to  use  near-horizontal  boundary  forces,  such  as  those  pro- 
duced by  tiebacks. 

COMPARISON  WITH  OTHER  PROGRAMS 

This  topic  has  been  approached  in  two  ways  (Siegel,  1975a;  Boutrup, 
1977).   The  first   is  to  simplify  the  problem  in  terms  of  boundary  and 
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subsurface  conditions,  and  to  compare  values  of  minimum  FOS  generated  by 
searching  with  a  variety  of  available  programs.  Since  only  STABL 
searches  with  shapes  other  than  circles,  the  comparison  comes  down  to 
values  generated  by  the  various  circular  methods  of  slices.  The 
Geotechnical  Engineering  literature  contains  many  such  comparisons, 
which  show  that  the  simplified  Bishop  assumption  used  in  CIRCL2  compares 
well  with  the  longer  and  more  iterative  total  equilibrium  approaches. 
Research  at  Purdue  shows  the  same  results. 

A  second  approach  is  to  generate  critical  surfaces,  in  relatively 
complex  problems,  by  using  the  non-circular  searching  modes  in  STABL. 
This  is  followed  by  generating  FOS  values,  for  the  specific  surface 
defined  by  STABL,  by  other  available  methods.  Putting  aside  the  proba- 
bility that  the  engineer  would  fail  to  find  the  most  critical  surface  by 
specifying  possibility  after  possibility,  the  simplified  Janbu  approach 
of  STABL  gives  somewhat  conservative  (low)  values  of  FOS. 

It  is  well  to  remember  two  facts  relative  to  comparisons.  Firstly, 
there  is  no  analysis  method  which  gives  the  correct  value  for  FOS  for 
reasonably  complex  problems.  Since  there  is  no  method  for  determining 
the  true  answer,  there  is  no  way  of  assigning  errors  to  any  practical 
method.  Secondly,  comparisons  vary  with  the  part  of  the  data  field 
being  examined,  i.e.,  relative  values  change  with  the  specific  set  of 
examples. 

Neither  of  the  two  difficulties  discussed  above  need  perplex  the 
engineer  using  STABL.   When  STABL  is  being  considered  as  a  replacement 
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for  extant  methods,  they  can  be  run  in  parallel  with  each  other,   and 
answers  compared. 

SUMMARY 

The  STABL  slope  stability  system  allows  computerized  searching  for 
a  variety  of  shapes  for  potential  sliding  surfaces.  All  such  surfaces 
are  represented  as  a  series  of  straight  lines,  and  the  line  directions 
are  chosen  in  a  pseudo-random  manner.  The  analyses  proper  use  simpli- 
fied Bishop  or  simplified  Janbu  methods  of  slices  to  determine  FOS 
values. 

Plotting  routines  are  necessary  to  assess  the  exact  shapes  and 
positions  of  generated  sliding  surfaces  and  to  determine  whether  it  is 
probable  that  the  more  critical  sliding  surfaces  have  been  generated  and 
analyzed.   Error  messages  allow  checks  on  input  and  execution  steps. 

The  STABL  system  should  be  entered  with  reasonable  estimates  of  the 
minimum  FOS,  so  that  iteration  of  the  implicit  expressions  for  FOS  is 
minimized.  The  estimates  may  be  obtained  from  existing  chart  solutions 
of  the  real  problem  in  simplified  form.  The  CIRCL2  searching  option  in 
STABL  should  yield  values  of  FOS  which  compare  closely  with  those  of 
other  circular  searches.  The  BLOCK  and  RANDOM  options  in  STABL  are 
likely  to  produce  somewhat  lower  values  of  FOS  than  specific  surface 
analyses  by  other  methods.  One  reason  for  this  difference  is  the  proba- 
bility that  a  limited  number  of  specific  surfaces  will  not  identify  the 
critical  sliding  surfaces.  The  STABL  searches  should  do  a  reasonable 
job  of  accomplishing  this  objective. 


21 


The  latest  documentation  relating  to  STABL  is  presented  in  Appendix 
B.  Appendices  E  and  F  contain  recommended  changes  in  the  program  with 
respect  to  the  earthquake  option,  and  the  convenient  entry  of  near- 
horizontal  tieback  forces  into  the  calculations  of  FOS. 
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APPENDIX  A 
EXAMPLE  PROBLEM 

The  stability  of  the  slope  shown  in  Figure  Al  was  investigated 
using  STABL.  The  example  subsurface  profile  consists  of  3  soil  types 
for  the  irregular  surface  and  circular  surface  analysis.  For  the  block 
analysis,  a  "weak"  seam  was  introduced,  as  soil  type  A,  to  illustrate 
this  option.  A  high  groundwater  level  was  assigned  to  the  slope. 

The  irregular  surface  search  was  confined  to  zones  where  the  ran- 
domly generated  surface  would  initiate  and  terminate  between  x  =  25  to 
45  ft  (7.62  to  13.72  m)  and  x  =  140  to  185  ft  (42.67  to  56.39  m), 
respectively.  The  50  surfaces  which  were  successfully  generated  between 
these  user-defined  limits  are  shown  in  Figure  A2.  A  summary  of  the  ten 
most  critical  surfaces  is  shown  in  Figure  A3  with  the  surface  represent- 
ing the  minimum  factor  of  safety  being  "highlighted"  with  asterisks  (*). 

A  similar  search  was  also  performed  for  50  circular  surfaces  with 
the  same  initiation  and  termination  limits  as  used  for  the  irregular 
surface  search.  The  50  surfaces  are  shown  in  Figure  A4  and  the  summary 
in  Figure  A5. 

The  "weak"  seam,  which  was  introduced  into  the  subsurface  profile, 
usually  leads  to  failure  represented  by  block  or  wedge  shapes.  The  cen- 
tral portion  of  the  block  failure  was  limited  between  the  two  boxes, 
judiciously  placed,  to  "force"  generation  of  block  surfaces  along  the 
"weak"  layer.   As  can  be  seen  from  Figure  A6,  the  major  portion  of  all 
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stratum  consisting  of  Soil  Type  h   vas  used  for  the  '"block'  analysis,  only. 
Figure  Al  Example  Slope  Geometry  and  Subsurface  Profile 

Table  Al  Summary  of  Soil  Properties  Used  in  Example 


Soil  Type  ft 

Total  Unit 
Wt.  (pcf) 

Sat.  Unit 
Wt.  (pcf) 

Cohesion 

Intercept 

(nsf) 

■ 
Friction 
Angle 
(deg) 

1 

125.0 

129.0 

0.0 

38.0 

2 

120.0 

127.0 

1+00.0 

21.0 

3 

128.0 

132.0 

600.0 

26.0 

h 

12  h .  0 

128.0 

200.0 

16.0 

For  S.I.  conversion,  use  : 

1  psf  =  1*7.873  Pascals  (or  1+ .  882  kg/sq.  metre) 
1  pcf  =  16.018  kg/  cuMc  metre 
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block  surfaces  was  confined  to  the  "weak"  layer.  Figure  A7   presents   a 
summary  of  the  most  critical  surfaces. 

The  three  most  critical  irregular,  circular  and  block  surfaces  were 
also  analyzed  further  for  earthquake  loads.  Seismic  coefficients  with  a 
value  of  0.02  were  used  with  the  vertical  loading  being  downwards  and 
the  horizontal  loading  outwards  from  the  slope.  The  change  in  factors 
of  safety  are  shown  in  Figures  A8  to  A10. 

Thus,  based  on  a  first  run,  we  can  summarize  the  minimum  factors  of 
safety  for  different  critical  surface  types: 


TYPE  OF 
CRITICAL  SURFACE 


FACTOR  OF  SAFETY 


STATIC  CASE 


EARTHQUAKE  LOADS 


Irregular 

Circular 

Block 


1.40 
1.54 
1.12 


1.30 
1.42 
1.03 


After  one  has  examined  the  output  of  the  initial  run,  the  user  can 
evaluate  the  extent  of  the  search  for  the  "most"  critical  surface  and 
decide  if  further  investigations  are  required. 
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APPENDIX  B 
STABL  DOCUMENTATION 

1.  STABL  Program,  available  on  unlabled  9-track  tape  for  IBM  or  CDC 
Fortran  codes  at  the  required: 

1.  recording  density,  in  bpi 

2.  fixed  block  size  (URL=80) 

3.  conversion  format  (EBCDIC  or  ASCII) 

2.  Listing  of  program 

3.  STABL  User  Manual,  JHRP-75-9 

4.  Computer  Analysis  of  General  Slope  Stability  Problems,  JHRP-75-8 

5.  Computerized  Slope  Stability  Analysis  for  Indiana  Highways,  Vol.  1, 
JHRP-77-25 

6.  Vol.  2  (1977  Program  Listing),  JHRP-77-26 

7.  Three-Dimensional  Slope  Stability  Analysis,  JHRP-81-17 

For  ordering  or  further  information,  please  contact: 

Sunil  Sharma  or  Prof.  C.  W.  Lovell 

Grissom  Hall,  School  of  Civil  Engineering 

Purdue  University 

West  Lafayette,  Indiana   47907 

Phone:   (317)  494-5025  (or  5034) 
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APPENDIX  C 
USERS  UNDER  FHWA  IMPLEMENTATION  PLAN 


Arkansas  State  Highways  and  Transportation  Department 

Materials  and  Research  Division 

P.  0.  Box  2261 

Little  Rock,  Arkansas  72203 

Attn  :  Mr.  Jake  Clements,  Jr. 

Phone:  501-569-2496 


2.   Colorado  Division  of  Highways 
4201  East  Arkansas  Avenue 
Denver,  Colorado  80222 
Attn  :  John  Gilmore 
Phone:  303-757-9275 


3.   Florida  Department  of  Transportation 
Soils  Materials  and  Research  Engineer 
Bureau  of  Materials  and  Research 
P.  0.  Box  1029 
Gainesville,  Florida  32602 
Attn  :  Dr.  H.  K.  Ho 
Phone:  904-372-5304 


4.  Georgia  Department  of  Transportation 
Data  Processing  Unit 
Materials  and  Research  Lab. 
Georgia  DOT 
15  Kennedy  Drive 
Forest  Park,  Georgia  30050 
Attn  :  Roger  Pruitt 
Phone:  404-363-7567 


Kansas  State  Department  of  Transportation 

Bureau  of  Computer  Services 

10th  Floor,  State  Office  Building 

Topeka,  Kansas  66612 

Attn  :  John  Lebbert 

Phone:  913-296-3657 


6.   Idaho  Department  of  Transportation 
Information  and  Data  Processing 
P.  0.  Box  7129 
Boise,  Idaho  83707 
Attn  :  Ron  Cole 
Phone:  208-334-2551 
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7.  Iowa  Department  of  Transportation 
Ames,  Iowa 

Attn  :  Bruce  L.  Brakke,  Division  Bridge  Engineer 
Phone: 


Louisiana  Department  of  Transportation 

P.  0.  Box  44245 

Baton  Rouge,  Louisiana   70804 

Attn  :  J.  B.  Esnard,  Soils  Design  Engineer 

Phone:  504-342-7623 


9.  Minnesota  Department  of  Transportation 
State  Transportation  Building 
St.  Paul,  Minnesota  55155 

Attn  :  Virgil  V.  Mikkelsen,  Foundation  Engineer 
Phone:  612-296-2304 


10.  Mississippi  State  Highway  Department 
P.  0.  Box  1850 
Jackson,  Mississippi  39205 
Attn  :  Mr.  L.  Ray  Marler 
Phone:  601-354-7211 


11.   Nebraska  Department  of  Roads 
P.  0.  Box  94759 
Lincoln,  Nebraska   68509 
Attn:   A.  R.  Kennedy 
Phone:   402-473-4752 


12.  Nevada  Department  of  Transportation 
Materials  and  Testing  Division 
1263  South  Stewart  Street 
Carson  City,  Nevada  89712 
Attn  :  Ted  Beeston 
Phone: 


13.  New  Jersey  Department  of  Transportation 
1035  Parkway  Avenue  CN600 
Trenton,  New  Jersey  08625 
Attn  :  Anil  Mehta 
Phone:  609-292-3574 
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14.   State  of  New  Mexico  Highway  Department 
Bureau  of  Materials 
P.  0.  Box  1149 

Santa  Fe,  New  Mexico  87504-1149 
Attn  :  Richard  Lueck 
Phone:  505-983-1149 


15.   North  Carolina  Department  of  Transportation 
Division  of  Highways 
P.  0.  Box  25201 

Raleigh,  North  Carolina  27611 
Attn  :  Mr.  Larry  Colbert 
Phone:  919-733-2075 


16.   Oklahoma  Department  of  Transportation 
Design  Support  Unit 
200  N.  E.  21st  Street,  Room  2C9 
Oklahoma  City,  Oklahoma   73105 
Attn  :  Charles  Whittle 
Phone : 


17.   Oregon  Department  of  Transportation 

Bridge  Section 

Oregon  DOT  Building 

Salem,  Oregon  97310 

Attn  :  John  Marks 

Phone:  503-378-6551 


18.  Texas  State  Department  of  Highways  and  Public  Transportation 
11th  and  Brazos 
Austin,  Texas   78701 
Attn  :  Mr.  Wayne  Henneberger 
Phone : 


19.  Turner-Fair bank  Highway  Research  Center 
6300  Georgetown  Pike 
McLean,  Virginia  22101 
Attn  :  Chien-Tan  Chang 
Phone:  703-285-2357 


20.   Utah  Department  of  Highways 
4501  South  2700  West 
Salt  Lake  City,  Utah  84119 
Attn  :  Loren  Rausher 
Phone:  801-965-4196 
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21.   Washington  State  Department  of  Transportation 
Materials  Laboratory 
P.  0.  Box  167 

Olympia,  Washington  98504 
Attn  :  Leroy  Wilson 
Phone:  206-753-7100 


22.   Wisconsin  Department  of  Transportation 
3502  Kinsman  Boulevard 
P.  0.  Box  7878 

Madison,  Wisconsin  53704-7878 
Attn  :  Gary  Whited 
Phone:  608-246-3249 


23.   Wyoming  Highway  Department 
Cheyenne,  Wyoming  82002 

Attn  :  Bill  Sherman,  Chief  Engineering  Geologist 
Phone:  307-777-7801 
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APPENDIX  D 
COMMENTS  ABOUT  SAFETY  FACTORS  FOR  SAMPLE  PROBLEMS 

Some  concern  has  been  expressed  about  the  discrepancy  between  the 
safety  factors  given  in  the  User  Manual  for  the  sample  problems  and 
those  computed  by  the  Users.  These  differences  are  due  to  use  of  STABL 
on  non-CDC  computers  and  thus,  we  would  like  to  offer  an  explanation  for 
these  variations. 

The  factors  of  safety,  for  surfaces  generated  using  random  numbers, 
will  depend  somewhat  on  the  sequence  of  such  numbers  and  the  method  of 
computation.  The  example  problems  presented  in  the  User  Manual  were 
analyzed  on  a  CDC  computer,  using  a  library  function  to  generate  the 
required  random  numbers.  The  library  function  available  for  VAX  comput- 
ers will  provide  the  same  sequence  of  random  numbers.  Thus,  the  results 
shown  in  the  example  problems  may  be  duplicated  only  if  the  program  is 
executed  on  a  CDC  or  VAX  computer.  For  other  cases  (e.g.,  IBM  comput- 
ers), the  function,  RANF,  within  the  program  source  code  is  used  to  gen- 
erate the  sequence  of  random  numbers.  However,  since  such  a  sequence  of 
random  numbers  will  be  different  from  the  one  generated  by  a  CDC  or  VAX 
library  function,  different  potential  failure  surfaces  will  be  analyzed 
by  the  program.  Consequently,  a  different  minimum  factor  of  safety  will 
be  computed  and  the  critical  surface  will  differ  somewhat  from  the  one 
presented  in  the  examples. 
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APPENDIX  E 
MODIFICATIONS  AND  REVISIONS  OF  STABL 

Some  additional  changes  have  been  made  in  the  last  18  months,  due 
to  a  much  more  extensive  use  of  STABL  for  teaching  purposes  and  also  as 
a  result  of  interaction  between  the  Users  and  Purdue  University.  Most 
of  these  changes  are  minor  and  only  improve  the  operation  of  the  pro- 
gram. However,  a  significant  modification  is  proposed  for  the  portion 
of  the  program  dealing  with  the  pseudo-static  earthquake  analysis. 

The  existing  program  lowered  factors  of  safety  (FOS)  for  the  effec- 
tive stress  analysis  where  excess,  positive  pore  pressures  were  con- 
sidered in  the  analysis.  As  mentioned  on  p.  14  of  the  User  Manual 
(JHRP-75-9),  (Ref.  10),  the  program  considers  the  generation  of  excess 
pore  pressures  based  on  a  seismic  coefficient.  These  pressures  are  com- 
puted according  to  the  inertial  forces  generated  by  the  seismic  coeffi- 
cient (see  p.  126,  JHRP-77-25,  Ref.  1).  We  feel  that  this  is  not  par- 
ticularly realistic,  especially  as  it  applies  only  to  effective  stress 
analysis.  Conventionally,  the  selected  seismic  coefficient  will  be 
expected  to  account  for  such  effects,  as  well  as  a  loss  of  strength  dur- 
ing an  earthquake. 

Thus,  we  have  revised  the  program  to  include  only  the  effects  of 
additional  inertial  forces  induced  by  the  specified  seismic  coeffi- 
cients, without  changing  the  static  pore  pressures.  For  this  revised 
method,  the  following  changes  are  proposed: 
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1.  Change   the   statement   on    line    WGHT    704    to    read: 
IF(RU(SOILTP).EQ.0..AND.CU(S0ILTP).EQ.0.)GO  TO  23 

2.  LINES  WGHT  708,  710,  712  should  be  deleted. 

Also  note  that  a  value  for  the  cavitation  pressure  is  required  when 
using  the  EQUAKE  option  (see  p.  41  in  User  Manual,  Ref.  10).  However, 
the  following  errata  should  be  inserted  in  the  User  Manual  to  replace 
the  second  paragraph  on  p.  14: 

The  inertial  forces  due  to  the  seismic  coefficients 
are  at  the  center  of  gravity  of  each  slice.  These  forces 
do  not  change  the  pre-earthquake  static  pore  pressures 
in  the  slope.   If  significant  excess  pore  pressures  changes 
or  loss  of  shear  strength  is  expected,  or  in  the  case  of  a 
"high-risk."  slope,  a  complete  dynamic  analysis  should  be 
performed. 

There  are  other  minor  changes  which  do  not  affect  the  FOS  computa- 
tions.  These  are  listed  below: 

1.  In  subroutine  PROFIL: 

Line  PR0F1151  should  read:   CALL  READER(DUMMY,NP(K),0) 

(date  of  change,  4-9-84) 

2.  In  subroutine  RANDOM: 

Line  RAND1130  should  read:   WRITE(6,110)ERROR(20) 

(date  of  change,  4-9-84) 

after  line  RAND1830,  insert  these  3  lines: 

IF(MB  .HE.  1)  WRITE(6,136) 
136   FORMAT(10X,'  *  *  SAFETY  FACTORS  ARE  CALCULTED  BY  THE  ', 
1    'MODIFIED  JANBU  METHOD  *  *',//) 

(date  of  change,  1-9-84) 
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3.  In  subroutine  RANSUF: 

Line  RANS  766  should  read:  IF(ILIMIT.EQ.O)  GO  TO  2 

Line  RANS  798  should  read:  GO  TO  8 

After  line  RANS1376,  insert:   IF(SURF(K,1) .LT.BPT)GO  TO  3 

After  line  RANS1464,  insert:   IF(SURF(K,1) .LT.BPT)GO  TO  3 

(date  of  change,  3-1-83) 

4.  In  subroutine  FACTR: 

After  line  FCTR  492,  insert:   IF(A1(I) .LT.0.0)A1(I)=0.0 

(date  of  change,  3-1-83) 

5.  For  those  Users  with  random  generator  subroutine  RANF, 

Line  RANF   4  should  read:   DATA  XN, SEED/131072 ,27487/ 

After  line  RANF   8,  insert:   IF(SEED.LT.O)SEED=SEED+XN 

(date  of  change,  3-1-83) 

NOTE:  All  of  the  above  are  FORTRAN  statements  and  should  be  inserted 
according  to  conventional  programming  rules  into  their  correct 
column  fields. 
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APPENDIX  F 
INPUT  OF  TIEBACK  LOADS 

F.l  Introduction 

The  use  of  tiebacks  in  geotechnical  enginering  and  construction  for 
stability  of  slopes  and  support  of  excavations  has  increased  substan- 
tially within  the  last  several  years.  As  a  result,  the  need  for  a 
method  of  analyzing  the  overall  stability  of  slopes  and  retaining  walls 
subjected  to  horizontal  or  inclined  concentrated  loads  has  become  more 
evident.  Until  now,  the  input  of  horizontal  or  inclined  concentrated 
loads  acting  on  a  near  vertical  slope  was  somewhat  difficult  in  STABL. 
In  addition  the  factor  of  safety  was  not  formulated  for  this  type  of 
loading  and  thus,  did  not  fully  account  for  the  distribution  of  force  to 
the  failure  surface  caused  by  concentrated  boundary  loads. 

Therefore,  to  increase  the  versatility  of  STABL,  new  routines  have 
been  created  within  STABL  to  permit  input  of  horizontal  or  inclined  con- 
centrated loads.  These  routines  were  created  specifically  for  the  input 
of  tieback  loads  but  may  be  easily  used  for  any  type  of  concentrated 
load  applied  to  the  ground  surface.  The  latest  version  of  STABL, 
STABL4,  contains  the  new  routines  which  utilize  Flamant's  Formulas  as 
proposed  by  Morlier  and  Tenier  (1982)  and  the  modified  Bishop  method  of 
analysis  for  circular  failure  surfaces,  and  the  simplified  Janbu  method 
of  analysis  for  non-circular  failure  surfaces.  The  tieback  option  may 
be  used  with  either  random  or  specific  failure  surface  generation 
methods  for  irregular,  block  or  circular  failure  surfaces.  Throughout 
this  appendix  and  within  STABL4  the  word  "tieback"  is  used  to  mean  tie- 
back  or  other  types  of  concentrated  loads  applied  to  the  ground  surface. 
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Tieback  or  other  types  of  concentrated  loads  are  input  by  specify- 
ing the  ground  surface  boundary  number  where  the  load  is  to  be  applied, 
the  Y  coordinate  of  the  point  of  application  on  the  ground  surface,  the 
magnitude  of  the  point  load,  the  horizontal  spacing  between  point  loads, 
the  inclination  of  the  load  as  measured  clockwise  from  the  horizontal 
axis,  and  the  length  of  the  tieback,  Figure  Fl .  For  concentrated  boun- 
dary loads  such  as  strut  loads  in  a  braced  excavation  which  do  not 
extend  into  the  ground  like  tiebacks,  the  length  of  the  tieback  is  zero. 
An  equivalent  line  load  is  calculated  for  each  tieback  load  specified 
assuming  a  uniform  distribution  of  load  horizontally  between  point 
loads. 

A  short  description  of  the  new  tieback  routines  is  presented  to 
help  the  User  understand  the  method  and  assumptions  used  in  STABLA  for 
analyzing  slopes  subjected  to  concentrated  loads.  Two  example  problems 
are  also  presented  to  demonstrate  the  input  of  horizontal  and  inclined 
loads  and  the  effect  of  concentrated  loads  on  the  factor  of  safety 
against  slope  failure.  The  example  problems  include  input  and  output 
data  and  plots.  In  addition,  the  input  format  for  concentrated  loads  is 
presented  along  with  input  restrictions,  error  codes,  and  a  listing  of 
the  additions  and  modificatons  to  STABL3  required  to  update  that  version 
to  STABL4. 

F.2  Description  of  New  Tieback  Routines 

Unlike  other  slope  stability  programs,  STABL4  distributes  the  force 
from  a  concentrated  load  throughout  the  soil  mass  to  the  whole  failure 
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surface  and  hence  to  all  slices  of  the  sliding  mass.  Other  slope  sta- 
bility programs  on  the  other  hand,  only  take  a  concentrated  load  into 
account  on  the  slice  on  which  it  acts.  This  distribution  of  load 
throughout  the  soil  mass  is  a  unique  feature  of  STABL4 . 

First  an  equivalent  line  load  is  calculated  for  a  row  of  tiebacks 
by  dividing  the  specified  tieback  load  (point  load)  by  the  corresponding 
horizontal  spacing  between  tieback  loads.  The  resulting  line  load  is 
called  TLOAD,  Figure  F2,  and  is  inclined  from  the  horizontal  by  an  angle 
INCLIN.  The  radial  stress  on  the  midpoint  of  a  slice  is  calculated 
using  Flamant's  Formula  (Morlier  and  Tenier,  1982): 


a  = 

r 


2(TL0AD)cos(TTHETA) 
(tO(DIST) 


where 


r 
TLOAD 

TTHETA 


DIST 


Radial  stress 

Equivalent  tieback  line  load 

Angle  between  the  line  of  action  of 
the  tieback  and  the  line  between  the  point 
of  application  of  the  tieback  on  the 
ground  surface  and  the  midpoint  of  a  slice. 


Pi 


Distance  between  the  point  of  application 
of  the  tieback  on  the  ground  surface  and  the 
midpoint  of  a  slice. 


The  radial  force,  PRAD,  at  the  midpoint  of  the  base  of  the  slice 
due  to  the  concentrated  load  is  calculated  by  multiplying  the  radial 
stress  by  the  length  of  the  base  of  the  slice: 


PRAD  => 


2 (TLOAD) cos (TTHETA) .    (DX) 


(it)  (DIST) 


cos (ALPHA) 
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where 

PRAD     =    Radial  force  on  base  of  slice 
due  to  concentrated  load 

ALPHA    =    Inclination  of  base  of  slice 

DX      =    Slice  width 

Note  that  the  radial  stress  produced  on  the  base  of  the  slice  by 
the  concentrated  load  is  proportional  to  the  load  applied  (TLOAD)  and 
the  width  of  the  slice  (DX),  inversely  proportional  to  the  distance 
between  the  point  of  application  of  the  load  and  the  midpoint  of  the 
base  of  the  slice  (DIST),  and  dependent  upon  the  angle  between  the  line 
of  action  of  the  load  and  the  line  between  the  point  of  application  of 
the  load  and  the  midpoint  of  the  base  of  the  slice  (TTHETA).  Therefore, 
slices  which  are  in  line  with  the  direction  of  the  concentrated  load 
will  receive  a  larger  portion  of  the  total  load  than  will  slices  which 
are  farther  away  and  whose  angle  TTHETA  is  large. 

The  radial  force  PRAD  is  distributed  in  the  same  manner  to  all  the 
slices  of  the  sliding  mass.  The  radial  forces  on  all  the  slices  are 
then  summed  in  the  direction  of  the  concentrated  load,  PSUM,  and  com- 
pared with  the  applied  load,  TLOAD.  Since  the  sum  of  radial  forces  for 
a  failure  surface,  PSUM,  is  not  always  exactly  equal  to  the  applied  load 
due  to  slope  geometry  and  the  shape  of  the  failure  surface,  the  radial 
force  applied  to  the  base  of  each  slice  is  modified  as  follows: 

PRAD  =  TLOAD/PSUM 

The  refined  radial  force  for  each  slice,  PRAD,  is  broken  into  its 
components  normal  and  tangential  to  the  base  of  the  slice  for  calcula- 
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tion  of  the  factor  of  safety.   The  normal  and  tangential  components   of 
the  force  due  to  the  concentrated  load  are  respectively: 

PNORM  =  (PRAD)cos(ALPHAl) 
PTAN  =  (PRAD)sin(ALPHAl) 

The  same  process  is  repeated  for  all  additional  rows  of  tiebacks. 
The  sum  of  the  normal  components  and  the  sum  of  the  tangential  com- 
ponents due  to  all  rows  of  tiebacks  are  then  used  in  the  slice  equili- 
brium equations  for  calculating  the  factor  of  safety. 

There  is  a  special  case  where  the  tieback  loads  will  not  be  distri- 
buted to  quite  all  the  slices  of  the  sliding  mass  and  is  shown  in  Figure 
F3.  Figure  F3  shows  the  limit  of  the  stress  distribution  for  a  benched 
slope.  The  force  due  to  the  applied  load  is  not  distributed  to  the 
slices  of  the  far  left  or  the  slices  of  the  far  right  since  this  would 
require  distribution  of  load  through  air  and  not  the  soil  mass. 

Examples  showing  the  use  of  the  new  tieback  routines  in  STABL4  are 
presented  in  the  next  section. 

F.3  Examples  Using  Concentrated  Loads 

The  details  of  the  examples  presented  are  not  of  concern  since  the 
primary  purpose  of  the  examples  is  to  demonstrate  the  use  of  the  new 
routines  for  concentrated  loads.  It  should  be  noted  that  the  critical 
failure  surface  and  the  minimum  factor  of  safety  will  change  in  most 
instances  due  to  the  introduction  of  concentrated  horizontal  or  inclined 
loads.  In  addition,  the  examples  presented  herein  only  analyze  a  speci- 
fied failure  surface  for  the  sake  of  simplicity  and  clarity.   A  random 
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search  should  always  be  performed  to  find  the  new  critical  failure  sur- 
face and  corresponding  minimum  factor  of  safety  against  overall  slope 
failure. 

The  profile  shown  in  Figure  F4  represents  a  proposed  excavation  in 
a  silty  sand  with  vertical  sides  supported  by  steel  sheet  piling.  It 
has  been  determined  that  the  critical  failure  surface  is  circular  and 
exits  the  slope  at  the  toe  of  the  sheeting.  The  minimum  factor  of 
safety  for  overall  stability  was  determined  to  be  0.451  for  the  critical 
surface  with  center  at  (41  ft,  76  ft)  (12.50  m,  23.16  m)  and  radius  41.9 
ft  (12.77  m). 

Two  methods  of  support  have  been  proposed.  The  first  method 
involves  using  two  rows  of  tiebacks  and  the  second  method  involves  using 
a  bracing  system  inside  the  excavation. 

F.3.1  Example  Problem  #1 

It  has  been  determined  that  two  rows  of  tiebacks  are  required  for 
local  stability  of  the  wall  (Figure  F5).  The  first  (upper)  row  of  tie- 
backs  will  be  horizontal  and  the  free  length  of  the  tiebacks  will  be  40 
ft  (12.19  m).  The  first  row  of  tiebacks  will  each  carry  a  load  of 
150,000  lbs  (68,040  kg)  and  will  be  spaced  10  ft  (3.05  m)  apart  horizon- 
tally. The  second  (lower)  row  of  tiebacks  will  be  inclined  15  degrees 
as  measured  clockwise  from  the  horizontal  and  will  have  a  free  length  of 
30  ft  (9.14  m).  Each  tieback  in  this  row  will  carry  a  load  of  200,000 
lbs  (90,720  kg)  and  will  also  be  at  a  10  ft  (3.05  m)  horizontal  spacing. 
Both  rows  of  tiebacks  will  transfer  their  loads  to  the  sheeting  using 
horizontal  steel  beams   for  wales.   Therefore,   the  assumption  of   a 
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uniform  distribution  of  load  between  the  tiebacks  will  be  reasonable  for 
this  case  since  the  action  of  the  wale  will  cause  an  equivalent  horizon- 
tal line  load  to  be  formed  for  each  row  of  tiebacks. 

It  is  desired  to  know  what  effect  the  tiebacks  have  on  the  factor 
of  safety  of  the  previously  mentioned  critical  failure  surface.  There- 
fore, the  program  will  be  run  using  a  specified  circular  failure  surface 
and  the  input  data  shown  in  Figure  F6.  The  output  for  this  run  is  given 
in  Figures  F7,  F8  and  F9. 

Note  that  if  the  unbonded  portion  of  any  tieback  does  not  intersect 
the  failue  surface,  i.e.,  terminates  within  the  sliding  mass  as  in  tie- 
backs  #1  and  #2  of  Figure  F10,  the  load  from  that  particular  row  of  tie- 
backs  will  not  be  considered  in  the  factor  of  safety  calculation.  Only 
the  load  from  tieback  #3  of  Figure  F10  would  be  considered  in  the  factor 
of  safety  calculation  since  the  failure  surface  intersects  the  unbonded 
portion  of  tieback  #3. 

As  will  be  seen  in  the  next  example,  the  length  of  a  "tieback"  for 
a  braced  system  of  support  is  equal  to  zero  and  all  concentrated  loads 
will  be  considered  for  all  failure  surfaces  in  the  factor  of  safety  cal- 
culation. 

F.3.2  Example  Problem  #2 

The  second  method  of  support  proposed  involves  using  bracing  to 
provide  local  and  overall  stability  to  the  wall  (Figure  Fll).  The  brac- 
ing system  will  consist  of  horizontal  struts  bearing  against  steel  wales 
which   transfer   their   loads   to   the  sheet  piling.   The  first  level  of 
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PROF I L 

TIEBACK  EXAMPLE  PROBLEM  *1 

4  4 

10.  35.  50,  35.  1 

50.  35.  50,1  60.  1 

50.1  60.  80.  62.  1 

80.  62.  93.  62,. 1 

SOIL 

1 

1.20,  125.  0.  18.  0,  0.  1 

WATER 

1  62.4 
4 

10.  35. 
50.  35, 
50.01  50. 
82.  51, 
TIES 
? 

2  53,  150000.  10,  0,  40, 
2  48.  200000.  10.  15,  30. 
SURBIS 

7 

50,  35,2 
56;  36,8 
64,  40,9 
70,  45.1 
74.  49. 
78,  54.1 
82.3  62. 
EXECUT 


FIGURE   F6  -    INPUT  DATA  FOR  EXAMPLE   //l 


—SLOPE  STABILITY  ANALYSIS — 
SIMPLIFIED  -IANPU  MFTHOD  OF  SI  ITES 
OR  SIMPLIFIED  BISHOP  METHOD 


PROBLEM  DESCRIPTION    TIEBACK  EXAMPLE  PROBLEM  *1 


BOUNDARY  COORDINATES 


4  TOP    BOUNDARIES 
4  TOTAL  BOUNDARIES 
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BOUNDARY 
NO. 


X-LEFT 
(FT) 

10.00 
50.00 
50.10 
80.00 


Y-LFFT 
(FT) 

35.00 
35.00 
60.00 
62.00 


X-RIGHT 
(FT) 

50.00 
50.10 
80,00 
93.00 


Y-RIGHT 
(FT) 

35.00 
60.00 
62.00 
62.00 


SOIL  TYPE 
BELOU  BND 

1 
1 
1 
1 


ISOTROPIC  SOIL  PARAMETERS 


1  TYPE(S)  OF  SOIL 


SOIL    TOTAL    SATURATED 

TYPE   UNIT  UT.   UNIT  UT. 

NO.    (PCF)      (PCF) 


1 


120.0 


1.25.0 


COHESION   FRICTION 


INTERCEPT 
(PSF) 


PORE     PRESSURE   PIEZOMETRIC 
SURFACE 
NO. 


ANGLE    PRESSURE   CONSTANT 
(DEG)   PARAMETER     (PSF) 


18.0 


0 


1 


1  PIEZOMETRIC  SURFACE(S)  HAVE  BEEN  SPECIFIED 
UNIT  WEIGHT  OF  WATER  =   62.40 

PIEZOMETRIC  SURFACE  NO.   1  SPECIFIED  BY   4  COORDINATE  POINTS 


POINT 

X-UATER 

Y-UATER 

NO. 

(FT) 

(FT) 

1 

10.00 

35.00 

2 

50.00 

35.00 

3 

50.01 

50.00 

4 

62.00 

51.00 

TIEBACK  LOAD(S) 

2  TIEBACK  LOAD(S)  SPECIFIED 


TIEBACK 
NO. 

X-POS 
(FT) 

Y-POS 
(FT) 

LOAD 
(LBS) 

SPACING 
(FT) 

INCLINATION 
(DEG) 

LENGTH 
(FT) 

1 
2 

50.0? 
50  .  05 

53.00 
48.00 

150000.0 
200000.0 

10.0 
10.0 

0 
15.00 

40.0 
30.0 

FIGURE  F7  -  OUTPUT  FOR  EXAMPLE  #1 


TRIAL  FAILURE  SURFACE  SPECIFIED  DY   7  COORDINATE  ROINTS 


POINT 

X-SURF 

Y-SURF 

MO. 

(FT) 

(FT) 

1 

50.00 

35.00 

2 

56.00 

36.80 

3 

64.00 

40.90 

4 

70.00 

45.  JO 

5 

74.00 

49.00 

6 

7P.00 

54-10 

7 

82.30 

62.00 
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CIRCLE  CENTER  AT  X 


41.0  J  Y 


76.0   ANH  RADIUS*    41.9 


FACTOR  OF  SAFETY  FOR  THE  PRECEDING  SPECIFIED  SURFACE  =   1.870 

WARNING  -  FACTOR  OF  SAFETY  IS  CALCULATED  FY  THE  MODIFIED  BISHOP 

METHOD.  THIS  MFTHOTi  IS  UALID  ONLY  IF  THE  FAILURE  SURFACE 
APPROXIMATES  A  CIRCLE. 


i  I  I  I  I  -t-  i  i  I  i  i  -t-  I  i  i  i  i 


I  I  I  I  +  I  I  I  I  i  +  i  I  I 


FIGURE   F8  -   OUTPUT   FOR  EXAMPLE    //l    (Continued) 
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struts  will  carry  150,000  lbs/strut  (68,040  kg/strut)  and  will  be  spaced 

10   ft   (3.05  m)   apart.   The  second  level  of  struts  will  carry  200,000 

lbs/strut  (90,720  kg/strut)  and  will  be  spaced  10  ft  (3.05  m)  apart. 

Again  the  assumption  of  a  uniform  distribution  load  between  the  struts 
is  reasonable  due  to  wale  action. 

Again  the  program  is  run  using  the  input  data  shown  in  Figure  F12 
for  the  braced  excavation  and  a  specified  failure  surface  and  the  output 
is  given  in  Figures  F13,  F1A  and  F15.  As  mentioned  in  Example  #1,  the 
length  of  a  "tieback"  for  a  concentrated  load  such  as  the  load  produced 
by  a  row  of  struts  is  input  as  zero.  All  concentrated  loads  with  length 
zero  will  be  used  for  the  calculation  of  the  factor  of  safety  for  all 
failure  surfaces. 


F.A.I   Data  Input  Format  for  Tieback  Loads 
(if  specified) 


COMMAND  CARD 
DATA  CARD 
DATA  CARD 


TIES      Command  Code 

Integer   Number  of  tieback  loads 

Integer   Boundary  number  where  tieback 
load  is  applied 

Real      Y  coordinate  of  the  point  of 
application  of  tieback  load 
(ft)  or  (m) 

Real      Load  per  tieback 
(lbs)  or  (kg) 

Real      Horizontal  spacing  between 
tiebacks  (ft)  or  (m) 

Real      Inclination  of  tieback  load  as 
measured  clockwise  from  the 
horizontal  plane  (deg) 

Real      Free  length  of  tieback  (ft)  or  (m) 
(Equal  to  zero  if  other  than  a 
tieback  load) 


Note:   Repeat  preceding  data  card  for  each  tieback  load. 
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PROF I L 

TIEBACK  EXAMPLE  PROBLEM  412 

4  4 

10,  35,  50,  35.  1 

50.  35.  50.1  60,  1 

50.1  60.  80.  62,  1 

80,  62.  93,  62,  1 

SOIL 


1 

120.  125.  0. 

18,  0 ,  0 . 

WATER 

1  62.4 

4 

10,  35. 

50,  35. 

50.01  50. 

82.  51. 

TIES 

2 

2  53,.  150000, 

10.0.  0 . 

2  48,  200000, 

10.  0.  0. 

SURBIS 

7 

50.  35.2 

56.  36.8 

64.  40.9 

70,  45,1 

74,  49, 

78.  54,1 

82.3  62. 

EXECUT 

FIGURE  F  12 -INPUT  DATA  FOR  EXAMPLE   #2 
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--'-:i  OPF  STARILTIY  ANALYSIS-- 
S.TMPI  IFTED  IANBU  METHOD  OF  ELITES 
OR  STMPL  IFTED  BISHOP  METHOD  " 


PROBLEM  DESCRIPTION    TIEBACK  EXAMPLE  PROBLEM  *? 


BOUNDARY  COORDINATES 

4  TOP    BOUNDARIES 
4  TOTAL  BOUNDARIES 

BOUNDARY      X-LEFT       Y-LEFT     X-R'IGHT     Y-RTGHT     SOU   TYPE 
NCI-  (FT)         (FT)         (FT)         (FT)       BELOW  DND 

1  10.00       35.00       50.00        35.00  1 

50.00        35.00       50.10        60.00  1 

3  50.10        60.00       80,00        62.00 


4  80.00        62.00        93.00        62.00  1 

ISOTROPIC  SOIL  PARAMETERS 

1  TYPE(S)  OF  SOIL 

SOIL    TOTAL    SATURATED    COHESION   FRICTION     PORE     PRESSURE   PIEZOMETRIC 
TYPE   UNIT  UT.    UNIT  WT.    INTERCEPT      ANGLE    PRFSSURE   CONSTANT      SURFACE 
NP-     <PCF>       (PCF)        (PSF)        (DEG)   PARAMETER      <PSF)  NO. 

1      120-0       125.0  0        18.0  0  O  ] 

1  PIEZOMETRIC  SURFACE(S)  HAUE  BEEN  SPECIFIED 
UNIT  WEIGHT  OF  WATER  =   62.40 

PIEZOMETRIC  SURFACE  NO.   ]  SPECIFIED  BY   4  COORDINATE  POINTS 


POINT 

X-UIATER 

Y-UATER 

NO. 

(FT) 

(FT) 

1 

10.00 

35,00 

2 

50.00 

35.00 

3 

50.01 

50.00 

4 

82.00 

51.00 

TIEBACK  LOAD(S) 

2  TIEBACK  LOAD(S)  SPECIFIED 


IEBACK 
NO. 

X-POS 

(FT) 

Y-POS 
(FT) 

LOAD 

(LBS) 

SPACING 
(FT) 

INCLINATION 
(DEG) 

LENGTH 
(FT) 

1 
2 

50.07 
50.  05 

53.00 
48.00 

150000.0 
200000.0 

10.0 
10.0 

0 
0 

0 

0 

FIGURE  F13  -  OUTPUT  FOR  EXAMPLE  #2 


rpIAL  FAILURE  SURFACE  SPECIFIED  BY   7  COORDINATE  POINTS 


p  n  t  f  j  T 

X-SURF 

Y-SURF 

un . 

(FT) 

(FT) 

1 

50.00 

35.00 

-i 

56.00 

36.80 

3 

64.00 

40.90 

4 

70.00 

45.  10 

5 

74.00 

49.00 

6 

78.00 

54.  10 

7 

82.30 

62.00 
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CIRCLE  CENTER  AT  X 


76.0   AND  RADIUS.    4  1.9 


FACTOR  OF  SAFETY  FOR  THE  PRECEDING  SPECIFIED  SURFACE  =   2.071 


DARNING  -  FACTOR  OF  SAFETY  IS  CALCULATED  BY  THE  MODIFIED  BISHOP 

METHOD.  THIS  METHOD  IS  VALID  ONLY  TF  THF  FAILURE  SURFACE 
APPROXIMATES  A  CIRCLE. 


-*-  I  1  I  1   I  +  I  i  I 


I   I  I  I   I  ■+■  I   I  I  1   I  ■+■  I   I  1  I   I  +  I   I  I  I   I  +  I   I  I  I  I  ■+■  I  I  I  I  I 


FIGURE  F1A  -  OUTPUT  FOR  EXAMPLE  #2  (Continued) 
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F.4.2   Input  for  Suppressing  or  Reactivating  Tieback  Loads 
(if  specified) 

COMMAND  CARD   TIES      Command  Code 

DATA  CARD      Integer   Number  zero  (0) 

F.5   Input  Restrictions 

1.  The  point  of  application  of  a  tieback  on  the  ground  surface  may  not 
be  at  a  ground  surface  boundary  node.  Use  a  slight  vertical 
offset,  (i.e.  70.01  instead  of  70). 

2.  No  more  than  10  tieback  loads  can  be  specified  but  they  can  be  in 
any  order. 

3.  The  inclination  of  a  tieback  must  be  equal  to  or  greater  than  zero 
degrees  and  less  than  90  degrees  as  measured  clockwise  from  the 
horizontal. 

4.  The  horizontal  spacing  between  tiebacks  must  be  greater  than  or 
equal  to  1  ft  (or  1  metre  if  using  SI  units). 

5.  The  length  of  a  tieback  must  be  equal  to  or  greater  than  zero  ft. 
Zero  is  used  for  loads  other  than  tieback  type  of  loads. 

F.6   TIES  Error  Codes 

TI01  -  An  attempt  has  been  made  to  suppress  or  reactivate  undefined 
tieback  loads.  Data  must  be  defined  by  a  prior  use  of  command 
TIES  before  they  can  be  suppressed.  Suppressed  data  can  not 
be   reactivated   if   command   PROFIL   has   been  used   in   the 
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execution  sequence  subsequent  to  their  use,  whether   the   data 

are  active  or  suppressed. 


TI02  -  The  number  of  tieback  loads  specified  exceeds  10.  The  prob- 
lem must  either  be  redefind  so  fewer  tieback  loads  are  used, 
or  dimensioning  of  the  program  must  be  increased  to  accomodate 
the  problem  as  defined. 

TI03  -  A  negative  coordinate  has  been  specified  for  the  tieback 
load  indicated  or  the  calculated  Y  coordinate  of  the  end  of 
the  tieback  is  negative.  All  problem  geometry  must  be  located 
within  the  1st  quadrant. 

TI04  -  The  inclination  limits  have  been  exceeded  for  the  tieback 
load  indicated.  The  inclination  of  a  tieback  load  must  be 
equal  to  or  greater  than  zero  (deg)  and  less  than  90  (deg)  as 
measured  clockwise  from  the  horizontal. 

TI05  -  The  point  of  application  of  the  tieback  load  specified  does 
not  lie  on  the  ground  surface  boundary  specified.  Check  the 
boundary  number  specified  and  the  Y  coordinate  of  the  point  of 
application  of  the  tieback  load  indicated. 

TI06  -  The  horizontal  spacing  between  tiebacks  for  the  row  of  tie- 
backs  indicated  is  incorrect.  The  horizontal  spacing  between 
tiebacks  must  be  greater  than  or  equal  to  1  ft  (or  1  metre  if 
using  SI  units). 

TI07  -  The  length  of  the  tieback  indicated  is  incorrect.  The 
length  of  a  tieback  must  be  greater  than  or  equal  to  zero 
(ft).  Zero  is  used  for  loads  other  than  tieback  type  of  loads. 
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F. 7   List  of  Program  Additions  and  Modifications 

The  following  pages  list  the  program  lines  which  have  been  modified 
and  new  lines  which  are  required  for  calculation  of  the  factor  of  safety 
considering  the  presence  of  concentrated  loads.  There  are  two  new  sub- 
routines, TIES  and  TRANS.  Program  lines  which  have  been  modified  for 
the  input  of  concentrated  loads  are  denoted  by  JRCMAR84  which  indicates 
that  they  were  modified  by  James  R.  Carpenter,  March  1984. 

The  last  two  pages  of  this  section  list  the  lines  which  must  be 
interchanged  in  and/or  added  to  STABL4  to  run  the  program  on  an  IBM  com- 
puter. 
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